Interchain interactions and magnetic properties of Li2Cu02 
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An effective Hamiltonian is constructed for an insuiating cuprate witli edge-sliaring cfiains Li2Cu02. 
Tlie Hamiltonian contains tlie nearest and next-nearest neigfiboring intracliain and zigzag-type in- 
terchain interactions. Tlie values of the interactions are obtained from the analysis of the magnetic 
susceptibility, and this system is found to be described as coupled frustrated chains. We calculate 
the dynamical spin correlation function ^(q, to) by using the exact diagonalization method, and show 
that the spectra of ^(q, u>) are characterized by the zigzag-type interchain interactions. The results 
of the recent inelastic neutron-scattering experiment are discussed in the light of the calculated 
spectra. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40. Gb 
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One-dimensional cuprates have received much atten- 
tion as reference systems of high-Tc superconductors 
with two-dimensional Cu02 planes. Recently, a variety 
of the compounds with edge-sharing chains, where CUO4 
tetragons are coupled by their edges, were synthesized 
and found to show unique physical properties. Li2Cu02 
is one of the typical compounds having such chains. As 
shown in Fig. 1(a), the chains run parallel to b-axis and 
are stacked along a- and c-axes.Q 

An important feature of the edge-sharing chains is that 
a nearest neighboring (NN) magnetic interaction Ji be- 
tween Cu spins strongly depends on Cu-O-Cu bond an- 
gle 6. In the case that 9 = 90°, the superexchange 
process via O ions, which contributes to antiferromag- 
netic (AFM) interaction, is suppressed due to the or- 
thogonality of Cu3d and 02p orbitals, and ferromag- 
netic (FM) contribution caused by, for example, direct 
exchange mechanispji- between Cu3c? and 02p orbitals, 
becomes dominant .Eru With increasing 9, the AFM su- 
perexchange interaction increases, and consequently Ji 
changes from FM to AFM interaction at a critical an- 
gle 9c. In the previous studyja 9c was estimated to be 
about 95° from the cluster calculation. For Li2Cu02 with 
9 — 94°, Ji was evaluated to be FM (<0) with magnitude 
of 100 K. In addition to Ji, a next-nearest neighboring 
(NNN) magnetic interaction J2, which comes from Cu-0- 
0-Cu path, also plays an important role in the magnetic 
properties. The interaction J2 is AFM (>0), and its mag- 
nitude is known to be comparable to | Ji|. Therefore, an 
appropriate model describing the edge-sharing chain is 
a spin 1/2 Heisenberg model with NN and NNN inter- 
actions (a J1-J2 model). The ground state of the Ji- J2 
model has been extensively studieds For J2/I Ji| < 1/4, 
it is a FM state, while for J2/I Ji| > 1/4, it is a frustrated 
state with incommensurate spin correlation. 

In Li2Cu02, AFM long-range order occurs at Tn—9 K, 
and the magnetic structure below -2^ is FM along a- 
and b-axes and AFM along c-axis.El The recent inelas- 
tic neutron-scattering experiment showed the existance 
of interchain (IC) interactions which bring about the 
ordering.El The analysis of the dispersions along a- and 



c-axes in the linear spin- wave theory revealed that the IC 
interaction between NN Cu spins is of the order of 10 K. 
The band calculation also showed the existence of large 
effective hoppings between neighboring chains, and thus 
the superexchange interactions.^ Following these facts, 
IC interaction plays an important role in the magnetic 
properties of Li2Cu02. The importance of IC interac- 
tions have also been pointed out in other edge-sharing 
cuprates, CuGeOs (Ref. 7) and Sri4Cu2404i (Refs. 8 
and 9). 

In this paper, paying attention to the IC interaction, 
we examine the magnetic properties of Li2Cu02 such as 
magnetic susceptibility and magnetic excitation by ap- 
plying the exact diagonalization method on finite size 
clusters. We construct a minimal model, taking into ac- 
count the IC interaction to the J1-J2 model. A set of 
parameter values is obtained from the analysis of the tem- 
perature dependence of the magnetic susceptibility x{T). 
We find that the compound is described as a system with 
frustrated chains coupled by zigzag-type IC interactions. 
In order to examine the magnetic excitation, we calcu- 
late the dynamical spin correlation function S'(q, uj) in 
the chain direction. The calculated spectra show a flat 
dispersion caused by the frustration due to J2 in the 
low energy region. On the other hand, in the high en- 
ergy region, there is a dispersion, the energy position of 
which corresponds to that obtained from the linear spin- 
wave theory. The spectra of the dispersion are, however, 
broad. We show that this dispersion is brought about 
by the IC interaction with zigzag-type structure. The 
experimental results are discussed in the light of our the- 
oretical results. 

We first construct a minimal model for Li2Cu02 in- 
cluding IC interaction. As shown in Fig. 1(a), Li ions are 
located between the chains. The IC interaction works in 
a- and c- directions via Li ions. Since the hatched chains 
in Fig. 1(a) are situated in b-c plane, the orbitals rele- 
vant to the electronic states in each chain are Gu3dy2_^2 
and 02py,z ones. The possible paths which give IC in- 
teractions along c-axis Jc and along a-axis Jq are shown 
in Fig. 1(a) by thick solid and dotted lines, respectively. 
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For Jc (solid line), the orbitals of Li ions (Lils, 2s) couple 
to 02px in one chain, but to Oipy z orbitals in another 
one. On the other hand, for Ja, they couple to only 02px 
orbitals in both chains. As a result, the magnitude of J a 
is expected to be much smaller than that of Jc due to or- 
thogonality of Cu3dy2_^2 and Qlp^ orbitals. In fact, the 
neutron-scattering experimental indicates that the width 
of magnon dispersion along a-axis is narrower than that 
along c-axis, and Ja is ~ 4 K, which is less than half of 
Jc- In the present study, we neglect J a for simplicity. 
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FIG. 1. (a) The crystal structure of Li2Cu02. The solid, 
open and hatched circles denote Cu, O and Li ions, respec- 
tively. The edge-sharing chains run along b-axis, and are 
stacked along a- and c-axes. The thin solid and open lines 
connect intra- and inter-chain neighboring ions, respectively. 
The thick solid and dotted lines represent paths which give 
the IC interactions along c- and a-axes via Li ions, respec- 
tively, (b) The schematic picture of Cu-O-Li-O-Cu paths. 
NN Cu ions [Cu(l) and Cu(2)] and NNN Cu ions [Cu(l) and 
Cu(3)] in c-direction are connected by the paths depicted by 
two solid and two dotted lines, respectively, (c) The Ji- J2- Jc 
model. The Cu sites with S=l/2 spin are depicted by the solid 
circles. Ji and J2 are interactions between NN and NNN Cu 
spins in chains. Two chains are coupled by IC interactions 
JcNN and JcNNN , which work between NN and NNN Cu spins 
in different chains, respectively. We take Jcnn= JcNnn=>/c. 
b = he^ and c = cec are the primitive vectors. The rectangle 
represented by the broken lines denotes a unit cell. 



For Jc, there are two paths connecting between NN 
Cu ions [Cu(l) and Cu(2)] and between NNN Cu ions 



[Cu(l) and Cu(3)] as shown by two solid and two dot- 
ted lines in Fig. 1(b), respectively. They give rise to 
zigzag-type IC interactions between NN Cu spins and 
between NNN ones, Jcnn and Jcnnn, respectively. From 
the consideration of the crystal structure and the orbital 
configuration, the magnitudes of JcNN and Jcnnn are 
found to be the same because each path makes an equal 
contribution to the two interactions. Therefore, we take 

^cNN = >/cNNN = Jc- 
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FIG. 2. The magnetic susceptibility x{T) of Li2Cu02. The 
square symbols denote the experimental results in the mag- 
netic field applied along the chains from Ref. 3. The solid 
curve is the theoretical result in the J\-J2-Jc model with 
Ji=— 100 K, J2=40 K and Jc=16 K, obtained numerically 
on the system with 8x2 sites, g-factor is taken to be 2 from 
Ref. 10. 



Based on the above consideration, we adopt a JTnJ2- 
Jc model shown in Fig. 1(c). In the previous study,H Ji 
was evaluated to be —100 K. We determine J2 and Jc 
based on the analysis of the experimental data of x{T)- 
By diagonalizing the Hamiltonian of the Ji- J2- Jc model 
for finite-size clusters with MxN sites (M is the num- 
ber of sites in a chain and N is the number of chains), 
the theoretical x(T) is obtained. We preliminarily cal- 
culated it in 4x2 and 4x4 clusters in order to check the 
size effect along the IC direction (N dependence), and 
found that the difference of xC^) between N=2 and 4 is 
so small that clusters with N=2 are enough to investi- 
gate the effect of the IC interaction on x(T). In Fig. 2, 
we show the result for a 8x 2 cluster by solid line together 
with the experimental data denoted by square symbols. 
A good agreement between theory and experiment is ob- 
tained for J2=40 K and Jc=16 K. This value of Jc is not 
far from the. value determined by the neutron-scattering 
experiments The theoretical results of x(^) reproduce 
well the divergent behavior of the experimental ones at 
low temperatures. We examined x(2^) of single chain 
(the J1-J2 model) in the previous study,u and obtained 
a peak in x{T) at T=40 K. Therefore, we conclude that 
x(r) is not explained unless Jc is not taken into account. 

Next, we calculate the dynamical spin correlation func- 
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tion ^(q, oj) in order to investigate the magnetic excita- 
tion. This is defined as follows. 

S{q,Lu) = ^ |(m|5^|0)pj(^ - + Eo), (1) 

m 

where |0) and |m) are the ground and excited states with 
energies Eq and Em, respectively. is the Fourier trans- 
form of spatial spin density. 




FIG. 3. (a) S{q,Lj) along the chain direction qb with 
Ji = -100 K, J2=40 K and Jc=0 K. The chain has 20 sites, 
(b) The same as (a) but Jc=16 K. The system has 12x2 sites, 
simulating the coupled edge-sharing chains in Li2Cu02. The 
S functions are convoluted with a Lorentzian broadening of 
0.3 meV. The thin solid curve denotes a magnon dispersion 
[Eq. (3)] obtained by the linear-spin-wave theory. 

It is instructive to consider S{q, us) for the single chain 
( Jc=0 case) in order to understand the effect of IC inter- 
action on the magnetic excitation. The calculation is per- 
formed by applying the exact diagonalization technique 
on a 20-site single chain with Ji=— 100 K and J2=40 K. 
Because the ratio of J2/I Ji|(=0.4) is more than 1/4, the 
system is described as a frustrated chain. The calculated 
results of S{q,uj) are shown in Fig. 3(a), in which qi, 
denotes a momentum of the chain direction. The spec- 
tra have a period of 27r, and are symmetric with respect 
to qb = TT. Since the ground state has incommensurate 



spin correlation, the spectral intensity is much larger at 
(7fc~l/27r and 3/27r than that at other momenta. Much 
of the spectral weight is concentrated in the very low en- 
ergy region (cj ^ 2 meV), and small weight spreads over 
the high energy region. 

We turn our attention to the system with finite Jc- 
Since we are interested in the effect of the IC interac- 
tion on the spectra of frustrated chain, S'(q, lo) along 
the chain direction is examined, ^(q, w)'s for 12x2 and 
6x4 clusters, where q is the two-dimensional vector with 
q = {qb,qc), are calculated to understand the effect of 
the IC interaction. By comparing results between the 
two clusters, we found that there is no essential differ- 
ence along the chain direction. Therefore, we consider 
S{(i,uj) for the 12x2 cluster. 

The operator for the 12x2 cluster is written as 

5^ = -i= ^ e'^^y^ (51., + 5|,,,e'(* + *)), (2) 

where S\ ^ and ^ are the z component of spin oper- 
ators at i-th. sites on A and B chains in Fig. 1(c), re- 
spectively. Hi is the y component of the position vector 
at z-th sites, and is the number of Cu sites. Since in 
the case of finite Jc, the system has two Cu sites in unit 
cell as shown in Fig. 1(c), the lattice period in the chain 
direction is changed from |b| to |b|/2. Correspondingly, 
the period of S'(q, us) is changed from 27r to 47r. 

The result of S{q,uj) for the 12x2 cluster is shown 
inJ'ig. 3(b), where the momentum of qc is taken to be 
O.llil The spectrum is completely different from that for 
the single chain shown in Fig. 3(a). The strong inten- 
sity is seen at qb ^ 2tt, and thus the spectrum becomes 
asymmetric with respect to qb = tt. The asymmetry is 
understood as follows. For example, consider the spectra 
at q = (0,0) and (27r,0). at these momenta is given 

^y^J:^iS%,, + Sh,^) and respec- 
tively. From these expressions, it is clear that the weight 
at q = (27r, 0) is larger than that at q = (0, 0), when the 
A and B chains are coupled antiferromagnetically. A flat 
dispersion with small intensity is seen around qb tt at 
^ 5 meV. This is a remnant of the continuum seen in 
Fig. 3(a). In addition, a dispersive spectrum with energy- 
maximum at qb — TT is seen in Fig. 3(b). This dispersion 
is consistent with that obtained by the linear spin-wave 
theory in which the intrachain FM and interchain AFM 
ordering is assumed; 

uj{qb) = ^/[Ji{cosqb - 1) + J2{cos2qb - 1) + SJJ^ 

— 64 J2 cos^ qb cos^ (3) 

Here, we note that Jc of zigzag-type is responsible for 
this dispersion. Jc connects four spins in a chain with 
one spin in the neighboring chain as shown in Fig. 1(c). 
Therefore, these four spins tend to align parallel so as 
to reduce the magnetic energy. This is why the IC in- 
teractions have large contribution to the dispersion. The 
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spectra of the dispersion are rather broad. This is clearly 
different from the single-peak structure in a simple FM 
chain. This is because FM alignment in the chain is dis- 
turbed by the quantum fluctuation caused by Jc and frus- 
tration by J2. 

Finally, we discuss the recent inelastic neutron- 
scattering data along b-axis for Li2Cu02 in the light of 
our theoretical results. The experiment has shown that 
(i) when the momentum is far from magnetic zone center 
qh = 27r, the spectral intensity is very small, and (ii) the 
lowest excitation appears at the magnetic zone center, 
and the dipersion has a minimum at gf, = vr.tl The fea- 
ture (i) is in good agreement with the momentum depen- 
dence of the spectral intensity shown in Fig. 3(b). t2l For 
the feature (ii), the structure at tt in the experiment may 
correspond to that aLw < 5 meV which is a remnant of 
the frustrated state.E^ On the other hand, the dispersion 
at w ^ 20 meV has not been observed experimentally. 
In order to find the dispersion, it will be necessary to do 
more detail experiment, especially, in the higher energy 
region. 

We have studied the magnetic excitation spectra of 
Li2Cu02, and found that the zigzag- type IC interac- 
tion brings about the dramatic difference seen in the 
spectra between Fig. 3(a) and (b). Therefore, it is 
meaningful to examine the spectra of other edge-sharing 
compounds such as La6Ca8Cu2404i and Ca2Y2Cu50io 
because they have different IC interactions, reflecting 
the difference of the crystal structures. For example, 
in La6Ca8Cu2404i, which has the structure similar to 
Sri4Cu2404i, IC interactions are expected to be smaller 
than that in Li2Cu02.EfEl We thus suppose that the 
weight of the dispersion at high energy region (^20 meV) 
is very small and much of the weight is concentrated in 
the low-energy region (^5 meV) in La6Ca8Cu2404i in 
contrast to Li2Cu02. A comparison among these mate- 
rials will yield a clearer understanding of the effect of IC 
interactions on the magnetic excitation. 

In summary, we have investigated the magnetic prop- 
erties of Li2Cu02 with edge-sharing chains, taking into 
account the interchain interaction. We determined the 
magnetic interactions in a chain and between chains by 
the analysis of the magnetic susceptibility. We also cal- 
culated the dynamical spin correlation function. The re- 
sults show a dispersion with broad spectra induced by 
the interchain interaction. The zigzag-type interchain 
interaction causes the dramatic effects on the magnetic 
properties of Li2Cu02. It is highly desirable that the in- 
elastic neutron-scattering experiment is performed in the 
wide energy region to understand the characteristics of 
the dispersion along the chain. 
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